Astronomy & Astrophysics manuscript no. Gd202 


February 2, 2008 


(DOI: will be inserted by hand later) 





o 
o 



> 
(N 

in 
in 
o 
^1- 
o 

Oh 
6 



Probing the Canis Major stellar over-density as due to the 

Galactic warp 

Y. Momany 1 , S. R. Zaggia 2 , P. Bonifacio 2 , G. Piotto 1 , F. De Angeli 1 , L. R. Bedin 1 , G. Carraro 1 

1 Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 2, 1-35122 Padova, Italy 
e-mail: momany , piotto , carraro , deangeli , bedin@pd . astro .it 

2 INAF - Osservatrio Astronomico di Trieste, via Tiepolo 11, 1-34131 Trieste, Italy 
e-mail: zaggia , bonifacioOts . astro . it 

Received April 21, 2004; accepted May 25, 2004 

Abstract. Proper-motion, star counts and photometric catalog simulations are used to explain the detected stellar over-density 
in the region of Canis Major, claimed to be the core of a disrupted dwarf galaxy (CMa, Martin et al. 2004 Bellazzini et 
al. I2004> . as due to the Galactic warp and flare in the external disk. We compare the kinematics of CMa M-giant selected 
sample with surrounding Galactic disk stars in the UCAC2 catalog and find no peculiar proper motion signature: CMa stars 
mimic thick disk kinematics. Moreover, when taking into account the Galactic warp and flare of the disk, 2MASS star count 
profiles reproduce the CMa stellar over-density. This star count analysis is confirmed by direct comparison with synthetic color- 
magnitude diagrams simulated with the Besancon models (Robin et al. 2003 I that include the warp and flare of the disk. The 
presented evidence casts doubt on the identification of the CMa over-density as the core of a disrupted Milky Way satellite. 
This however does not make clear the origin of over-densities responsible for the ring structure in the anticenter direction of the 
Galactic halo (Newberg et al. 2002 , Yanny et al. 2003; Zaggia et al. 2004, in preparation). 
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1. Introduction 



Recent large scale surveys, in the optical and near-infrared, 
have been excellent tools to constrain the structure and the star 
formation history of the Milky Way and its satellite system. 
Growing evidence of satellite accretion, stellar streams, and 
sub-structures all point to inhomogeneities in both the Galactic 
disk and halo. The Sgr dwarf spheroidal (Ibata et al. I1997I is 
a strong evidence of the hierarchical formation in galaxies like 
the Milky Way. Ever since its discovery, the search for extra- 
Galactic satellite remnants has been most appealing. With the 
availability of the Sloan Digital Sky Survey, Newberg et al. 
( 2002 1 showed evidence for at least 5 over-densities in the anti- 
center direction of the Galaxy. Four of these features, possibly 
part of a ring, were found close to the Galactic plane, suggest- 
ing two possible origins: (1) the remnant of a dwarf satellite 
galaxy in the process of disruption; or (2) a particular distri- 
bution of Galactic disk stars. The ring structure has been con- 
firmed kinematically by Yanny et al. (2003 1 and using 2MASS 
data by Rocha-Pinto et al. (2003 1. Frinchaboy et al. (2004 1 and 
Crane et al. ( 2003 1 showed the existence of a structure of stellar 
clusters with coordinated radial velocities, further confirming 
the presence of a possible stellar ring at a distance of ^ 18 kpc 
in the anti-center direction. 

Martin et al. (|2004 M04), investigating the ring structure 
with the 2MASS all-sky catalog, pointed to an elliptical-shaped 



stellar over-density centered at / = 240°, b = -7°. Accounting 
for simulations of (a) past mergers thick disk formation (Abadi 
et al. (2003 1, and (b) dwarf galaxies in-plane accretion (Helmi 
et al. 2003 1, M04 proposed this "feature" as the core of a satel- 
lite galaxy currently undergoing an in-plane accretion, namely 
the Canis Major dwarf spheroidal galaxy (CMa). In a compan- 
ion paper, the same group (Bellazzini et al. 2004 1 searched for 
photometric signatures of CMa in color-magnitude diagrams 
(CMDs), and best identified the red clump and red giant branch 
of CMa in the CMDs of the Galactic open clusters NGC 2477 
and Tom 1. Forbes, Strader & Brodie (2004i studying the age- 
metallicity relation for CMa probable cluster members pointed 
out a clearly distinct extragalactic origin of CMa debris. On the 
other hand, Kinman, Saha & Pier (2004i did not find any over- 
density of RR Lyrae stars in an anti-center field of the Galaxy. 
This contrasts with the Zinn et al. d2"0TJ3 1 study who possibly 
identified CMa stars in the so-called southern arc of the CMa 
ring. 

Accreting Galactic satellites could be responsible for the 
creation of the observed Milky Way warp (Castro-Rodriguez 
et al. 2002 1. Evidence for the existence of a Galactic warp for 
the inter-stellar gas and dust are long dated since early HI ob- 
servations (Oort, Kerr & Westerhout 1958 1, and more recent 
inter-stellar diffuse dust emission (Freudenreich et al. U994l and 
Freudenreich 1996 1 studies. The warp has been also detected 
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in stellar star counts and can be seen as a systematic variation 
of the mean disk latitude with longitude (Djorgovski & Sosin 
1989 1: observationally, the warp is an upward bending from the 
Galactic plane in the first and second galactic longitude quad- 
rants (0° < I < 180°), and downward in the third and fourth 
ones (180° < / < 360°). The warp feature is therefore impor- 
tant in the direction of CMa and has to be accounted for when 
asymmetries across the Galactic plane are investigated. 

A comprehensive formation scenario for the features found 
in CMa and their connection with other galactic structures se- 
riously call for the determination of basic observational prop- 
erties of CMa such as: metallicity, proper motion (pm), radial 
velocity etc. In the present paper we address the mean pm of 
stars in the direction of CMa in search for any kinematic sig- 
nature. We also investigate a possible connection of the CMa 
over-density with the warp and flare of the Galaxy (Robin et 
al. 2003 and references therein) which manifest their maximum 
extent in the CMa region. 

2. CMa proper-motion 

In the absence of specific pm studies of CMa, in this paper we 
make use of the recent USNO CCD Astrograph Catalog second 
release (UCAC2, Zacharias et al. 2004 1 to investigate a possible 
pm signature of the CMa. The first epoch data for the pm come 
from different catalogs including Hipparcos/Tycho, AC2002.2, 
as well as re-measured AGK2, NPM and SPM plates. In gen- 
eral, Zacharias et al. estimate pm errors to be 1 4- 3 mas/yr for 
\2 ,h magnitude stars and about 4-^7 mas/yr around the 16' /! 
magnitude. The pm reference is absolute and it has been eval- 
uated using external galaxies. A useful feature of the UCAC2 
catalog is the inclusion of the J, H and K magnitudes from the 
2MASS survey. We also made use of the on-line Milky Way 
photometric and kinematic simulations based on the Besancon 
Galaxy model (Robin et al. 2003 1: its main advantage is the re- 
duction of a number of free parameters (e.g. the thin disk scale 
height at different ages) and, most interestingly, the inclusion 
of warping and flaring of the Galactic disk. 

To estimate the pm of CMa we extract UCAC2 pa, pd (in 
this paper we assume pa = pacosd), J, and K data in 3 fields 
having a radius of 1°. The first centered on the Sgr to be used 
for testing the procedure, the second on CMa and the last on the 
open cluster NGC 2477. Adopting the selection criteria out- 
lined in Bellazzini et al. and M04, we use the K, [J - K) di- 
agrams to extract M-giants presumably belonging to CMa and 
estimate the mean UCAC2 pm. Although simple, this approach 
is rather tricky. As it is well depicted by the study of Platais et 
al. {2003 1, the presence of residual color/magnitude terms in 
any pm analysis is almost unavoidable. This is specially true 
when one deals with an all-sky catalog, where the pm is de- 
rived from different first epoch sources. Aware of these facts, 
the minimal approach was to (i) assess the reliability of the 
UCAC2 measurements by a comparison with fields of known 
pm, and (ii) limit the UCAC2 pm data to only one first epoch 
catalog in order to limit internal inhomogeneities. Several tests 
were performed re-measuring the mean pm of globular clusters 
and comparing these with the values reported in Dinescu et al. 
dl999> . The overall agreement has been very good for all mea- 
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Fig. 1. 2MASS K, ( J - K) dereddend CMDs and corresponding 
UCAC2 proper-motion plots of the Sgr dwarf spheroidal, CMa 
and NGC 2477. For each row, from left to right: (1) a 1° Z, 
(J - K) 2MASS CMD limited to stars with UCAC2 pm data 
based on the AC2002.2 catalog; (2) the same CMD but show- 
ing only the selected stellar populations for which a mean pm 
is to be determined; and (3) the mean pm of the selected popu- 
lations following a 2<x-clipping to eliminate outliers (see online 
colored version of the Figure). 

surable cases (16 clusters with distances ranging from 3 to 13 
kpc), with an overall scatter of less than 2.5 mas. In conclusion, 
the best confidence level was found with the AC2002.2 dataset 
(see Zacharias et al. 2004): we restricted our analysis on it. 

Fig. n presents CMDs and pm diagrams for 3 selected 
fields, from top to bottom: Sgr-center, CMa-center (I = 240°, 
b = -7°) and NGC2477 -field. The stellar populations selected 
in the middle panels are mainly field stars chosen to span all 
possible color and magnitude ranges in order to check for any 
pm gradient. The bulk of the selected disk populations with 
(J-K) < 0.8 are expected to have a heliocentric distance of 
< 2 kpc (open triangles), while objects with (J-K) > 0.8 
and 8 < K < 14 are found mainly at distances > 2 kpc (filled 
circles). On the top row of Fig.[T]the identification of the Sgr 
M-giants (filled circles) is straightforward, showing the red gi- 
ant tip at K si 11 .0. The right upper panel plots the estimated 
mean pm of the 6 presumably ^zeW populations (filled triangles) 
which show a clustering around < pa, pd >^ (-1.0, -6.0) mas. 
The dispersion of the mean pm of the field stars is due to the 
different nature and distance of the selected objects. The dis- 
persion of the 6 field samples (^ 2.0 mas) can be used as a 
conservative error estimate of our mean pm determinations. 
Relatively speaking, the mean pm of the Sgr red giants clearly 
stands out at < pa,p5 >^ (-2.5, -2.0) mas, and is well sepa- 
rated from surrounding field stars. This is due to the almost po- 
lar orbit of Sgr. A filled square with error bars shows the mean 
value of the Sgr proper motion as derived by Irwin et al. ( 1996 1 
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from HST data: < \m, fid >= (-2.80 + 0.80, - 1 .40 + 0.80) mas. 
The excellent agreement between the two values strengthens 
our confidence in measuring pm with UCAC2 up to a distance 
of ~ 25 kpc from the Sun. 

Turning our attention to CMa, we extract M-giants within 
an oblique stripe, as done in Bellazzini et (2004i in the NGC 
2477 field. CMa giants are plotted as filled circles in the cen- 
tral panels of the middle and lower rows of Fig. Stars be- 
longing to the open cluster were selected within a radius of 
0°2 (light starred symbols), clearly showing the main sequence 
and red clump. As in the Sgr field, the mean pm of disk stel- 
lar populations (filled triangles) in the CMa-center and NGC 
2477 fields show a dispersion of +2.0 mas, with a rather lower 
dispersion around the NGC 2477 as expected for its low lat- 
itude. We estimate a mean pm of CMa-center M-giants to be 
< fia,fiS >- (-1.7,3.1) mas. This is in the same pm direc- 
tion of disk stars, although slightly offsetted (~ 1.5 mas in both 
fia and /jS). To better understand this offset, on both right pan- 
els we overplot contour levels of simulated kinematic stellar 
catalogs using the Besancon Galactic model. The simulated 
catalogs were divided into two stellar subgroups (as done for 
the 2MASS J, K color selection) according to their heliocen- 
tric distance: thin line contours show stars within 2 kpc, i.e. 
(J - K) < 0.8, while thick line contours show stars with dis- 
tance > 2 kpc, i.e. (J-K) > 0.8. For each population, the counts 
have been divided in 9 linear contour levels. We note that the 
above stated value of CMa-center pm is in fair agreement with 
expected mean pm value of the latter group of disk objects, 
with a mean distance of 8 kpc. Thus, when accounting for 
a mean error of +2.0 mas and an offset merely due to distance 
projections, the mean pm of CMa-center seems perfectly com- 
patible with outer disk dynamics in circular prograde motion. 
This finding is in agreement with the Crane et al. (2003) so- 
lution for the kinematics of the "Monoceros" structure which 
appears compatible with a disk component rotating at 220 km/s 
at a distance of 1 8 kpc and showing a thick disk velocity disper- 
sion. Radial velocity measurements and a detailed computation 
of the CMa orbit will be presented in a separate paper (Zaggia 
et al., in preparation). 

3. CMa and the Galactic disk warp 

Ldpez-Corredoira et al. ( 2002 1 recently studied the existence of 
the warp in the Galactic stellar populations using star counts of 
old-population red clump giants selected from 2MASS CMDs. 
Their Fig. 15, on the ratio of North and South Galactic caps 
star counts shows a clear sinusoidal behavior with longitude. 
Excess of star counts in the northern hemisphere latitudes was 
found for / < 180°, and an opposite trend for I > 180°. They 
basically concluded their analysis confirming the existence of 
the warp in the old stellar population whose amplitude is co- 
incident with that of warped gas and young disk stars. They 
also identified a strong flare, i.e. a change in the scale-height, 
of the outer disk. Previously, Freudenreich et al. (19943 us- 
ing COBE/DIRBE near-infrared data found that the absolute 
maximum amplitude of the latitude displacement is 1°5 + 2°0. 
Interestingly, the over-density in CMa and "structure A " in the 
symmetric quadrant (0° < / < 180°) (see Figs. 4 and 5 in M04) 
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Fig. 2. Latitude profiles of a 2MASS M-giant star sample se- 
lected around CMa. Upper-left panel assumes a North/South 
symmetry around b — 0°, lower-left panel assumes a warp am- 
plitude of 2°0 in the southern direction. Star counts in left pan- 
els are corrected for reddening using the Schlegel et al. maps. 
Right panels show the same plots, except for correcting the 
Schlegel et al. values with the formula given in Bonifacio et 
al. 



coincide respectively with the southern and northern regions 
where the warp amplitude is strongest. This clearly indicates 
that the warp and flare of the disk can heavily affect star counts 
analysis of low latitude, distant regions. To test the impact of 
Galactic warp and flare near the core of CMa we revisited the 
analysis performed by M04. A suitable CMa M-giants sample 
was extracted from the 2MASS catalog using the same oblique 
selection region depicted in the CMDs of Fig. 1, i.e. sampling 
the region with CMa giants excess. The sample was spatially 
limited along a latitude stripe between -20° < b < 20°, in 
the longitude range 235° < I < 245°, i.e. corresponding to the 
CMa spatial FWHM, as derived in M04. J and K magnitudes 
were corrected for reddening in two different ways using: (1) 
raw Schlegel et al. ( 1998 1 values, and (2) modified Schlegel et 
al. values as suggested in Bonifacio, Monai & Beers (2000 1. 
Star counts were computed in latitude bins of 0°5 (the results 
don't change using the 4° x 2° binning of M04). 

In Fig. |2 the upper-left panel shows the folded star count 
profiles for the northern (dashed line) and southern (solid line) 
galactic caps. In this panel the assumed North/South symmetry 
is at b — 0° without introducing any warping in the Galaxy, i.e. 
with a warping angle w = 0.0. Also plotted are the Schlegel et 
al. (1998) reddening mean values for the southern (filled rect- 
angles) and northern (open rectangles) hemispheres along the 
same spatial region. The CMa location is identified with a ver- 
tical arrow. A glance at Fig. |2] shows that the southern part of 
the profile is clearly offsetted with respect to the northern part. 
In the lower-left panel of Fig.|2]the same star counts are pre- 
sented assuming a warp amplitude of w — 2°.Q in the southern 
direction, i.e. a North/South symmetry at b = -2°0. The curves 
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plotted in the lower panel show no significant changes varying 
w by +0°3 around w = — 2°0. The clearest feature in this panel 
is the almost perfect symmetry around b = -T.Q, Moreover, 
the North/South symmetry is also reflected in the two redden- 
ing curves which show a far better symmetry down to the mid 
plane. As a result of taking into account the disk warp, one 
notes that although the CMa over-density persists, it is seri- 
ously weakened. 

Various authors have noted that the Schlegel et al. maps 
overestimate reddening for E(B - V) > 0.2, in particular, 
Bonifacio et al. (2000 1 proposed a linear correction [their for- 
mula (1)] which lowers asintotically E(B - V) of 35%. 
The correction becomes quite important at low latitude with 
a difference of as much as =2 0.2 mag in (J - K) for an 
E(B - V) = 1 .0. The previous analysis is repeated applying the 
Bonifacio et al. correction (right panels in Fig.|2j>. Assuming a 
w = -2°0 and correcting the reddening values it is quite evi- 
dent that both star counts and reddening profiles show an ex- 
cellent symmetry. This time however, the CMa over-density has 
almost disappeared. The assumed value of w is in good agree- 
ment with that found by Freudenreich et al. ( 1994, their Fig. 3) 
and Freudenreich ( 1996) for this particular zone. These studies 
suggest that the layers of dust and neutral hydrogen are simi- 
larly displaced from the Galactic plane. As regarding to stellar 
warp, evidence of a similar w value has been suggested in the 
2MASS analysis of Lopez-Corredoira et al. (20021. Indeed, for 
the line of sight at I = 240° and galactocentric distances in the 
range 10 -f 14 kpc they predict the height of the mid plane (due 
to the warp) to range between 300 - 450 pc. This translates 
into a warping angle of w — L8 -f 2° 4 in the South direction. 
In conclusion, an appropriate consideration for the warp and 
flare completely erases any over-density in CMa. 

Further support of the importance of the warp comes from 
the comparison of observed and synthetic CMDs. In Fig.[3]we 
compare a V, (B-I) CMD of NGC 2477 (upper panel, Momany 
et al. 2001) with a Besancon synthetic catalog (middle panel). 
There is an excellent agreement between all observed features 
and simulated ones, except of course for the absence of the 
open cluster component and differential reddening effects in 
the simulated CMD. Clearly, sequences previously attributed 
to CMa (blue plume, young main sequence, red giants and red 
clump) are fully reproduced by the Galactic warped model. 
Plotted as heavy symbols are stars with a heliocentric distance 
between 7-^9 kpc, a mean metallicity of [Fe/H]= -0.45 ± 0.25 
and a mean age of 5.0 ± 1.5 Gyr. Curiously, similar parame- 
ters have been derived in Bellazzini et al. (2004) for the CMa 
stellar populations: a mean heliocentric distance of 8.3 ±1.2 
kpc, a mean metallicity in the range -0.7 <[M/H]< 0.0 and a 
mean age between 2 ~ 7 Gyr. On the other hand, the Galactic 
warped model shows a considerable decrease (a ~10 factor) 
of the stellar component between 7 -j- 9 kpc in the NGC2477- 
North field (Fig. 3, lower panel), i.e. it reproduces quite well the 
CMa South/North over-density. Similarly, Besancon simulated 
diagrams reproduce CMDs of the background populations of 
the open clusters having -18° < b < -14° (see Sect. 4.3 in 
Bellazzini et al. 2003), showing a substantial decrease of the 
warp background contribution. 




B-I 

Fig. 3. V, (B - I) diagrams of NGC2477. Upper panel is from 
Momany et al. d2001t . middle panel is a Besancon simulated 
CMD in an equivalent area along the same line of sight (b = 
-5°8). Lower panel is a Besancon simulated CMD at positive 
latitude (b = +5°8). 

In conclusion, accounting for the Galactic warp and flare 
explains the detection of both the CMa over-density and related 
stellar populations. As a consequence, the CMa feature can not 
be the progenitor of the ring found by Newberg et al. 
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